home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / zrotg.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  830 b   |  31 lines
  1.       SUBROUTINE ZROTG( CA, CB, C, S )
  2. *     .. Scalar Arguments ..
  3.       DOUBLE PRECISION  C
  4.       COMPLEX*16        CA, CB, S
  5. *     ..
  6. *     .. Local Scalars ..
  7.       DOUBLE PRECISION  NORM, SCALE
  8.       COMPLEX*16        ALPHA
  9. *     ..
  10. *     .. Intrinsic Functions ..
  11.       INTRINSIC         CDABS, DCMPLX, DCONJG, DSQRT
  12. *     ..
  13. *     .. Executable Statements ..
  14.       IF( CDABS( CA ).NE.0.0D0 )
  15.      $   GO TO 10
  16.       C = 0.0D0
  17.       S = ( 1.0D0, 0.0D0 )
  18.       CA = CB
  19.       GO TO 20
  20.    10 CONTINUE
  21.       SCALE = CDABS( CA ) + CDABS( CB )
  22.       NORM = SCALE*DSQRT( ( CDABS( CA / DCMPLX( SCALE, 0.0D0 ) ) )**2+
  23.      $       ( CDABS( CB / DCMPLX( SCALE, 0.0D0 ) ) )**2 )
  24.       ALPHA = CA / CDABS( CA )
  25.       C = CDABS( CA ) / NORM
  26.       S = ALPHA*DCONJG( CB ) / NORM
  27.       CA = ALPHA*NORM
  28.    20 CONTINUE
  29.       RETURN
  30.       END
  31.